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Abstract 

This paper presents the numerical simulations of confined three dimensional coaxial water jets. The 
objectives are to validate the newly proposed nonlinear turbulence models of momentum and scalar 
transport, and to evaluate the newly introduced scalar APDF and DWFDF equation along with its 
Eulerian implementation in the National Combustion Code (NCC). Simulations conducted include the 
steady RANS, the unsteady RANS (URANS), and the time-filtered Navier-Stokes (TFNS) with and 
without invoking the APDF or DWFDF equation. When the APDF or DWFDF equation is invoked, the 
simulations are of a hybrid nature, i.e., the transport equations of energy and species are replaced by the 
APDF or DWFDF equation. Results of simulations are compared with the available experimental data. 
Some positive impacts of the nonlinear turbulence models and the Eulerian scalar APDF and DWFDF 
approach are observed. 

1.0 Introduction 

In this study we have focused on two subjects. The first one is to validate the newly proposed 
nonlinear models of turbulent momentum and scalar transport implemented in the NCC code (Ref. 1) 
using the experimental data from confined swirling coaxial water jets (Ref. 2), and pay special attention to 
the relative performance of the nonlinear scalar flux model vs. the linear scalar flux model. The second 
one is to validate the newly introduced scalar APDF (Ref. 3) and DWFDF (Ref. 4) equations and their 
Eulerian implementations in the NCC code using the same experimental data. 

Both validations have been carried out with three types of numerical simulation approaches. The first 
two approaches are the steady RANS and the unsteady RANS (URANS). The third one is the time 
filtered Navier-Stokes (TFNS (Refs. 5 and 6) approach. 

The experimental study provided flow structures that are resembling those often existing in a gas turbine 
combustor, for example, center swirling recirculation near the front of combustor, massive swirled 
separations at the front comers, strong swirling flow extended all the way to the exit of combustor, and 
significant changes of concentration in the radial direction at the front of combustor, etc. Detailed 
experimental data on the velocity and scalar concentration distributions is available for validating the 
turbulence models and evaluating the numerical simulation approaches. Since water was used in the 
experiment but NCC solves the compressible Navier-Stokes equations supplemented by the equation of state 
for ideal gas, the water flow experiment was converted to its corresponding air flow simulation according to 
the Reynolds number similarity mle under the condition of low speed. The results of air flow simulations 
were then rescaled back to their water flow counterparts and compared with the experimental data. 

The turbulence models for momentum and scalar transport validated in the current simulations are 
listed in Section 2.1. The APDF and DWFDF equation invoked in the hybrid approach of RANS/APDF, 
URANS/APDF and TFNS/DWFDF simulations are listed in Section 2.2. 
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The results of validations of the turbulence models are presented in Sections 3.1, 3.2 and 3.3. The 
results obtained from the linear scalar flux model and those from the nonlinear scalar flux model in all 
three types of simulations (RANS, URANS and TFNS) are compared with the available experimental 
data. The results obtained by invoking the APDF and DWFDF equation are presented in a similar way in 
Sections 3.4, 3.5 and 3.6. These simulations are of hybrid nature, i.e., the velocity field is solved by the 
continuity and momentum equations, but the transport equations of energy and species are replaced by the 
equation of APDF or DWFDF, and the APDF or DWDFD equation is solved by an Eulerian Monte Carlo 
particle method (Ref. 7). Finally, a brief summary is given in the conclusion. 

In addition, an Appendix is presented to compare the different performance of using a standard k-& 
eddy viscosity linear stress model vs. using the present k- 8 nonlinear stress model in the RANS 
simulations under the exact same conditions. This comparison along with the experimental data clearly 
demonstrates the improvements made by the turbulent nonlinear stress model. 


2.0 Turbulence Models and Scalar PDF Equation 

The models for turbulent stress tensor Xy and scalar flux vector 0 Z as well as the scalar APDF and 
DWFDF equation employed in the current simulations are listed in this section. 


2.1 Linear Model Versus Nonlinear Model 

2.1.1 Linear Models of Turbulent Stresses and Scalar Fluxes 

Based on the Boussinesq approximation of the linear relationship between the turbulent stress and the 
strain rate of flow field, the linear model of turbulent stress is formulated as 


1 

*ij ~~ 8 ij x kk =- 2 / C n P— (-% -Sy%/3) 


Similarly, the linear model of turbulent scalar flux is formulated as 


©i =-flr 


0p0 

dxi 


( 2 . 1 ) 


( 2 . 2 ) 


2.1.2 Nonlinear Models of Turbulent Stresses and Scalar Fluxes 

A general constitutive relationship between the turbulent stresses Xy and the strain rate of flow field 
Sjj , Cly suggests (Refs. 6 and 8) 


T ij ~\h x *k = -2fC il p^-(s iJ -8j,-5 tt /3) 

-A /pC(s; v A y - n lk s kj ) (2.3) 

+2 A fp—^Cl ik S^j - sf k £l kj + Cl ik S km Q. m j - £l k iSi m £l mk 8 jj / 3 + II s (Sy - 8 yS kk /3) , 

where, Sy = {Uij +U j,i)/2, fly ={Utj -U =(s kk S mm -S kl S lk )j 2 . The model coefficients C^A 3 

and A 5 are constrained by the realizability condition and the rapid distortion theory. They are formulated 
as (see Ref. 9): 
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in which, 


4=V6coscp, (p = ^-arccos(V61f*), W* (2.5) 

3 v 7 (o / 

C/*=V(V) 2 +(Q*) 2 , s*=js^sfj, Q*=^VV, (2.6) 

The coefficient /is a function of filtering Resolution Control Parameter (RCP) that is defined as a ratio of 
the time filter width A 7 T 0 a global integral time scale of the flow T\ RCP = A t !T and 


Ar | [ 


/ ^ *2 — - — 
l T ) { T { T 


As discussed in Reference 6 , RCP may be viewed as a percentage measure of the unresolved subscale 
turbulent kinetic energy relative to the total turbulent kinetic energy. Therefore, the value of RCP and the 
coefficient / are always between 0 and 1, and /= 1.0 in RANS and URANS simulations,/ < 1.0 in TFNS 
simulations. 

Similarly, the nonlinear model for scalar fluxes is formulated as (Ref. 10), 


®i ~ ~&T ~( c i S v + c 2 Cly )- 


Where denotes the turbulent diffusivity for the corresponding scalar quantity 0. It is often 
approximated by = v/Pr 9 , where Pr e represents the turbulent Prandtl number or Schmidt number 
depending on whether the scalar quantity 0 is the internal energy e or the species mass fraction . 

The turbulent eddy viscosity is defined as v T =f ■ C M • £ 2 /s. The coefficients c x and c 2 in the current 
simulations are set to be c x = c 2 =-02A. 

2.1.3 Model Equations of Turbulent Kinetic Energy and Its Dissipation Rate 

The (subscale or total) turbulent kinetic energy and its dissipation rate k , 8 will be determined from 
the following model equations: 


d d_„ 

— p/v H P Hf k = 

dt dx; 


i k = 4~ (P + Pr)/^ -4% -P £ 

dxj \_ dxj 


^ ^ ^ 5 SO 8^ 

— P £ +— pM,-e = — (n + n r )— — e -C sl x ij s ij --C E2 K — (2.U 

dt oxf uXf |_ oxf \ k k 

where C 8i and C £2 are model coefficients. We have adopted the commonly used values of C e i = 1.45 and 
C e2 = 1.92 in the present simulations. 
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2.2 Scalar APDF and DWFDF Equation 


The transport equation for the scalar APDF or DWFDF, I\\ t (\|/; .v. i ) , can be written as (Refs. 3 and 4) 


dF 0 | d(UiF<b) 

St SXf 


{a 


\ dF <s>)\ 

[dxi 


dxj )\ 


+ /• 


d f 1 

dVk\ 


\ 

) 


k = l,2,---,M + \ 


( 2 . 11 ) 


where /is defined in Equation (2.7) and 


l 

x 




S ij S ij + Q/y Clij 


X >sJ(v + V T )/s 


( 2 . 12 ) 


3.0 Results of Numerical Simulations 

Figure 1 shows the computational domain and the unstructured grid (Ref. 11), which consists of 
849,189 tetrahedral elements having 152,555 nodes. Based on the Reynolds number similarity, the case of 
swirling coaxial water jets was scaled to the case of swirling coaxial air jets, so that it is appropriate to use 
the NCC code to simulate this low Mach number flow. The water dye concentration was represented by 
the mass fraction of the dyed air. There are two inlet boundaries, one is for the inner tube, at which the 
velocity, temperature and density were specified: 25.87 m/s, 300 K, and 1.1774 kg/m 3 ; another one is for 
the annulus, at which radial profile was specified. The annular air passes through eight swirlers resulting 
in an annular swirling jet dumped into a suddenly expanded chamber and mixed with the inner dyed air 
jet. Complex flow features exist in the chamber: shear layers, massive separations near the front corner of 
the chamber and the center recirculation, etc., see Figure 2. At the exit of the chamber, the pressure of 

1.0 atm. was specified for RANS and URANS simulations, whereas an unsteady convective boundary 
condition was specified for TFNS simulations. 

The results of numerical simulations will be presented in the following six sections, each one focuses 
on a different simulation approach, namely, RANS, URANS, TFNS, RANS/APDF, URANS/APDF, and 
TFNS/DWFDF. 
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Figure 1. — Computational domain and grid. 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
URANS-sp_L, 2nd=-.01, 4th=.05 
dt=1.e-05, 18,000 time step. 
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Figure 2. — Flow structures in the center X-Y plane. 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1 .0, 2nd=-.01, 4th=.05 
RANS at 1 20,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1 .0, 2nd=-.01 , 4th=.05 
RANS at 1 20,000 iteration 
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Figure 3. — Global feature of RANS with linear and nonlinear scalar flux models. 


In the first three sections, the main interest is in the turbulent scalar flux models, i.e., the performance 
of nonlinear scalar flux model versus the performance of linear scalar flux model. It is noted here that the 
turbulent stress model used in the current simulations is always the nonlinear model that has already been 
validated in the past studies. In this study we use the experimental data of dye concentration to validate 
the newly proposed nonlinear scalar flux model. 

In the next three sections, the main interest is in the hybrid approach, in which the scalar APDF or 
DWFDF equation (Refs. 3 and 4) is invoked to replace the transport equations of energy and species, and 
an Eulerian solver for the APDF or DWFDF equation is employed. The numerical results from the hybrid 
RANS/APDF, URANS/APDF and TFNS/DWFDF approaches are compared with the experimental data. 
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3.1 Results of RANS Simulations With Linear and Nonlinear Scalar Flux Models 

The parameter setting in NCC for the steady RANS calculations are: cfl = 1.0, mass-imbalance = 

1 .0x10 3 , 2nd-smoothing-coeff = -0.01, 4th-smoothing-coeff = 0.05. The global feature of the simulations 
is shown in Figure 3, the left one is RANS with the linear scalar flux model, the right one is RANS with 
the nonlinear scalar flux model. More detailed results and comparisons will follow. 

3.1.1 Evolutions of Velocity Components at Probes Along the Centerline 

Figure 4 shows the iteration history of velocity components at probes 3 and 5. They indicate that the 
computed flow field is not perfectly steady (see the axial velocity u ), especially near the front part of the 
chamber. 


RAN S -sp_L RAN S -sp_nL 




Figure 4. — Iteration history of the computed velocity components at probes 3 and 5. 
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RANS, axial velocity profile along centerline 



Figure 5. — Centerline distributions of axial velocity and concentration. 


3.1.2 Centerline Distribution of the Axial Velocity and the Dye Concentration 

Figure 5 indicates that results from both simulations (using linear as well as nonlinear turbulent scalar 
flux model, respectively) are consistent with the experimental data, and there are no appreciable 
differences between the results from the two models. 

3.1.3 Radial Distribution of the Concentration at Downstream Locations 

The radial profiles of dye concentration at five downstream locations (x = 13, 25, 51, 102, 203 mm) 
are shown in Figure 6. Again, results from both simulations match the experimental data reasonably well, 
except for the location at x = 5 1 mm. The differences between the results from the two models are 
negligible. 

3.1.4 Radial Distribution of Mean Velocity Components f7, V 9 W at Downstream Locations 

The radial profiles of axial velocity [/, radial velocity V and azimuthal velocity W at seven 
downstream locations (x = 5, 25, 51, 102, 152, 203, 406 mm) are presented in Figure 7, Figure 8 and 
Figure 9. Results from both simulations match the experimental data quite well, except for the azimuthal 
velocity W at the further downstream locations x = 203, 420 mm, where the computed W are under 
predicted compared with the experimental data. The differences of results from the two models are 
negligible. 
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Concentration Concentration 




Figure 6. — Radial profiles of mean concentration at downstream locations. 
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Figure 7. — Radial profiles of axial velocity U at downstream locations. 
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Figure 7. — Concluded. 


NAS A/TM— 20 14-218134 


10 


v (m/s) v (m/s) 






Figure 8. — Radial profiles of radial velocity V at downstream locations. 
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RANS, Radial mean velocity 
X = 203 mm, at 120,000 


RANS, Radial mean velocity 
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RANS, Azimuthal mean velocity 



RANS, Azimuthal mean velocity profile 
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RANS, Azimuthal mean velocity profile 





Figure 9. — Radial profiles of azimuthal velocity W at downstream locations. 
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RANS, Azimuthal mean velociry profile 
15 p X=406mm, at 120,000 iteration 
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Figure 9. — Concluded. 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 
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Figure 10. — Comparison of center recirculation zone. 


3.1.5 Flow Structure and Contour Plots of Variables in the X-Y Center Plane 

To further examine the effects of turbulent scalar flux models (linear vs. nonlinear) on the flow and 
concentration features, we have plotted the center recirculation zone visualized by the iso-surface of zero 
axial velocity (see Figure 10), the contour plots of velocity components (Figure 1 1), the dye mass fraction 
(Figure 12), the temperature, density, gauge pressure and Mach number (Figure 13), the turbulent kinetic 
energy and its dissipation rate, and the effective viscosity (Figure 14). The comparisons of these various 
quantities indicate that, for the current steady RANS simulations of nonreacting flow, the differences in 
effects due to the linear and nonlinear scalar flux model are negligible. 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens =0.01 
Sp_L, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 




3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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Figure 11. — Comparisons of U,V, W contours in center plane. 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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Figure 12. — Comparisons of O 2 mass fraction of the dye in center plane. 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 


Y 



T 

3 . 002E+02 
3 . 001E+02 
3. 001E+02 
3 . 001E+02 
3 . 00 OE-f 02 
2 . 999E+02 
2 . 999E+02 
2 . 999E+02 
2 . 998E+02 
2 . 998E+02 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1 .0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 


! 2 . 720E+02 
2 . 057E+02 
1.393E+02 
7 . 300E+01 
6 . 667E+00 
-5 . 967E+01 
-1.260E+02 
-1 . 923E+02 
-2 . 587E+02 
-3 .250E+02 



3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1 .0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 
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Figure 13. — Comparisons of temperature, density, pressure and Mach number in center plane. 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 
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Inner-dye-jet, Turb Intens = 0.01 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_L, cfl=1.0, 2nd=-.01 , 4th=.05 
RANS at 120,000 iteration 

I 7 . 695E-03 
5 . 776E-03 
3 . 857E-03 
■ 1 . 938E-03 
® 1.866E-05 


mu 

1 . 729E-02 
1 . 537E-02 
1 . 345E-02 
1 . 153E-02 
9 . 614E-03 



3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
Sp_nl, cfl=1.0, 2nd=-.01, 4th=.05 
RANS at 120,000 iteration 
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Figure 14. — Comparisons of turbulent kinetic energy, dissipation rate and effective viscosity in center plane. 


3.1.6 Remarks 

1 . Water coaxial jets are satisfactorily simulated by air coaxial jets based on the Reynolds number 
similarity and under the low speed (nearly incompressible) condition. 

2. Concentration of the dye (issued from the inner jet) is well calculated by using both of the linear and 
the nonlinear species flux models, except in a narrow region near the front stagnation point. 

3. Calculations of other flow quantities are also satisfactory: 

a. Axial and radial velocity components ( U and V) are close to the experimental data; however the 
azimuthal velocity decays too fast towards downstream and significantly lower than the 
experimental data in the rear half part of the flow. 

b. The center recirculation bubble is shorter than the measurement by 25 percent. 

4. The nonlinear species flux model has not shown appreciable effects on the concentration field 
comparing with the results from the linear model. 
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3.2 Results of URANS Simulations With Linear and Nonlinear Scalar Flux Models 


In this section, we compare the unsteady RANS (URANS) results obtained from the linear scalar flux 
model with the results from the nonlinear model. The URANS calculations use the RANS solutions as 
their initial conditions. The numerical parameter setting for the NCC are: dt=1.0xl0“ 5 , cfl=2.0, inner 
iteration convergence-criteria^l.OxlO 1 , 2nd-smoothing-coeff = -0.01, and 4th-smoothing-coeff = 0.05. 
When the simulation becomes statistically steady, the mean flow variables are established by averaging 
over the last 10,000 time steps of the calculation. 

The global feature of the simulations is shown in Figure 15, the left one is URANS with the linear 
scalar flux model, the right one is URANS with the nonlinear scalar flux model. More detailed results and 
comparisons follow. 


3.2.1 Centerline Distributions of Mean Axial Velocity and the Dye Concentration 

Figure 16 indicate that both models yield results which are consistent with the experimental data, and 
the differences of their results are not appreciable. 


3D Roback-Johnson flow 
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URANS-sp_L, 2nd=-.01 , 4th=. 05 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
URANS-sp_NL, 2nd=-.05, 4th=.09 
dt=1.e-05, 18,000 time step. 
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Figure 1 5. — Global features of URANS with linear vs. nonlinear scalar flux models. 




Figure 16. — Centerline distribution of mean axial velocity U and mean dye concentration. 
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3.2.2 Radial Profiles of the Mean Concentration at Downstream Locations 

The radial profiles of the dye concentration at six downstream locations (x = 13, 25, 51, 102, 152, and 
203 mm) are shown in Figure 17. Again, the results from both models match the experimental data 
reasonably well, except for at the location x = 5 1 mm, and the differences among their results are 
negligible. 







Figure 17. — Radial distributions of dye concentration at downstream locations. 
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3.2.3 Radial Distribution of Mean Velocity Components U 9 V, W at Downstream Locations 

The radial profiles of axial velocity £/, radial velocity V and azimuthal velocity W at eight 
downstream locations (x = 5, 25, 51, 102, 152, 203, 305, 406 mm) are presented in Figure 18, Figure 19 
and Figure 20. Results from both models match the experimental data quite well, except for the azimuthal 
velocity W at the further downstream locations x = 203 - 406 mm, where the computed W are under 
predicted compared with the experimental data. 





Figure 18. — Radial distributions of mean axial velocity at downstream locations. 


NAS A/TM— 20 14-218134 


20 





3 

2.5 
2 

1.5 
► 1 


Radial Profiles 
X = 406 mm 


O Exp. 44 

A Exp. 45 

SP_L 

SP_NL 


0.5 - 
0 - 


_aa@- — ©ta a"A'zX36»GL©0^ 


■© — &GQ 


-0.5 - 


_1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 

r/R 


Figure 18. — Concluded. 
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Figure 19. — Concluded. 
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Figure 20. — Comparison of azimuthal velocity component W at downstream locations. 
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Figure 20. — Concluded. 


3.2.4 Flow Feature and Contour Plots of Variables in the X-Y Center Plane 

To further examine the effects of turbulent scalar flux models on the flow features, we have plotted 
the center recirculation zone visualized by the iso-surface of zero axial velocity (Figure 21), the contour 
plots of velocity components (Figure 22), the dye mass fraction and turbulent kinetic energy (Figure 23), 
the density and temperature (Figure 24). The comparisons of these variables indicate that, for the current 
URANS simulations of nonreacting flow, the differences of results due to the two scalar flux models are 
negligible. 
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3D Roback-Johnson flow 
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URANS-sp_L, 2nd=-.01, 4th=.05 
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Figure 21 . — Comparison of center recirculation zone. 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
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3D Roback-Johnson flow 
Inner-dye-jet Turb Intens = 0.01 
URANS-sp_nL, 2nd=-.01, 4th=.05 
dt=1.e-05, 12,000 time step. 
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3D Roback-Johnson flow 
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Figure 22. — Comparison of mean velocities U, V, W contours in center plane. 
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Figure 23. — Comparison of mean O 2 mass fraction of the dye and k contours in center plane. 
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Figure 24. — Comparison of mean density and temperature contours in center plane. 
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3.2.5 Remarks 

1 . Concentration of the dye (issued from the inner jet) is satisfactorily calculated by both of the linear 
and nonlinear species flux models, except in a narrow region near the front stagnation point. 

2. Other flow quantities also are satisfactorily calculated: 

a. Axial and radial velocity components (U and V) are close to the experimental data; however the 
azimuthal velocity decays too fast towards downstream and significantly lower than the 
experimental data in the far downstream of the chamber. 

b. The center recirculation bubble is shorter than its experimental counter part by 25%. 

3. The results have not shown appreciable effect of the nonlinear species flux model on the 
concentration field when compared with the results from the linear model. 

3.3 Results of TFNS Simulations With Linear and Nonlinear Scalar Flux Models 

In this section, we present and compare the results of TFNS simulations using linear and nonlinear 
subscale scalar flux models. The TFNS calculations use the URANS solutions as the initial condition. The 
numerical parameter setting in the NCC is: dt = 2.0xl0“ 6 , cfl=1.0, inner iteration convergence-criteria = 
l.OxlO -1 , 2nd-smoothing-coeff = 0.0, and 4th-smoothing-coeff = 0.05. At the outlet, an unsteady 
convective boundary condition is employed. The filtering resolution control parameter RCP is set to be 
0.25. The same grid used in the RANS and URANS simulations is used for the current TFNS simulations. 

The global feature of the TFNS simulations are shown in Figure 25, the left one is TFNS with the 
linear subscale scalar flux model, the right one is TFNS with the nonlinear subscale scalar flux model. 
More detailed results and comparisons follow. 

3.3.1 Time Variation of Instantaneous Velocity Components at Centerline Probes 

Time histories of velocity components at five downstream locations are recorded during the 
simulations, the last three of probes are shown in Figure 26. The ones on the left are TFNS with the linear 
subscale scalar flux model; the ones on the right are TFNS with the nonlinear subscale scalar flux model. 
It can be seen that the turbulent fluctuation quickly develops from the location of probe 3 to the region of 
probes 4 and 5, so that the simulated flow in the chamber is a fully developed turbulent flow. 
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Figure 25. — Global features of TFNS using linear and nonlinear subscale scalar flux models. 
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Figure 26. — Time history of velocity components u, vand wat probes 3, 4 and 5. 
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3.3.2 Centerline Distributions of Mean Axial Velocity and Mean Dye Concentration 

Figure 27 indicate that results from both scalar flux models are consistent with the experimental data, 
and the differences between their results are small, albeit appreciable. 

3.3.3 Radial Profiles of Mean Concentration at Downstream Locations 

The radial profiles of mean concentration of the dye at six downstream locations (x = 13, 25, 51, 102, 
152, 203 mm) are shown in Figure 28. Again, results from both models match the experimental data quite 
well, including those at the location of x = 51 mm, and this is a noticeable improvement over the RANS 
as well as the URANS results. There are small but appreciable differences between the results from the 
two subscale scalar flux models. 



Figure 27. — Centerline distribution of mean axial velocity L/and mean concentration of the dye. 



Figure 28. — Radial distributions of mean dye concentration at downstream locations. 
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Figure 28. — Concluded. 


3.3.4 Radial Profiles of Mean Velocity Components at Downstream Locations 

The mean radial profiles of axial velocity U, radial velocity V and azimuthal velocity W at eight 
downstream locations (x = 5, 25, 51, 102, 152, 203, 305, and 406 mm) are presented in Figure 29, Figure 
30 and Figure 31. Results from both models match the experimental data quite well. For the azimuthal 
velocity W, results from the TFNS simulations show significant improvement over the results from the 
RANS and the URANS simulations near the wall region of the front chamber (x=5, 25,51 mm) as well as 
at the further downstream locations x = 203 - 406 mm. Again, the differences between the results from 
the two subscale scalar flux models are small. 
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Figure 29. — Radial distributions of mean axial velocity at downstream locations. 
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Figure 29. — Concluded. 



Figure 30. — Radial distributions of mean radial velocity at downstream locations. 
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Figure 30. — Concluded. 
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Figure 31 . — Radial distributions of mean azimuthal velocity at downstream locations. 
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Figure 31. — Concluded. 
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Figure 32. — Comparison of the center recirculation zone. 


3.3.5 Mean Flow Structure and Contour Plots of Variables in the X-Y Center Plane 

To further examine the effects of the subscale scalar flux models on the mean flow quantities, we 
have plotted the center recirculation zone visualized by the iso-surface of zero mean axial velocity (Figure 
32), the contour plots of mean velocity components and vorticity magnitude (Figure 33), the mean 0 2 
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mass fraction of the dye, the subscale kinetic energy and its dissipation rate (Figure 34). These 
comparisons indicate that results from these two models are very close. 
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Figure 33. — Comparison of mean velocities U, V, W and vorticity contours in center plane. 
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Figure 34. — Distributions of O 2 mass fraction of the dye, subscale kinetic energy and dissipation rate in 
center plane. 


3.3.6 Instantaneous Contour Plots of Variables in the X-Y Center Plane 

Additional comparisons in terms of the instantaneous features are presented here. Figure 35 shows the 
center recirculation zone visualized by the iso-surface of zero axial velocity, the contour plots of 
instantaneous velocity components, the vorticity magnitude, the 0 2 mass fraction of the dye, the subscale 
kinetic energy, and the effective viscosity p + p r . Comparison of these snap shots suggests that the 
nonlinear subscale scalar flux model adds finer details to the unsteady flow features. 
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Figure 35. — Contours of instantaneous variables in X-Y center plane. 
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Figure 35. — Concluded. 
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3.3.7 Remarks 


1. Concentration of the dye (issued from the inner jet) is well calculated at all streamwise locations. 

a. Other flow quantities are also satisfactorily calculated: All velocity components (U,V,W) were in 
good agreement with the experimental data, except the azimuthal velocity decays faster towards 
downstream when compared with the experimental data. The center recirculation bubble matches 
its experimental counterpart. 

b. Results from the TFNS simulations show significant improvements over the results from the 
RANS and URANS simulations. But they have not suggested a dramatic effect of the nonlinear 
subscale species flux model on either the concentration field or the flow field, i.e., the effect is 
positive, but small or very small. 


3.4 Results of RANS Simulation Invoking the Scalar APDF Equation 

In this section, the results of a steady RANS simulation invoking the scalar APDF equation are 
presented. The scalar APDF equation is solved by an Eulerian method. The main focus is to compare the 
performance of this hybrid approach with the performance of the pure RANS approach. 

Figure 36 shows the contour plots of the axial velocity and the 0 2 mass fraction of the dye in the X-Y 
center plane, and they are almost identical to their RANS counterparts presented in Section 3.1. Figure 37 
and Figure 38 present the comparisons of the results from the hybrid RANS/APDF simulation and from 
the pure RANS simulation at downstream locations x = 51 and 102 mm. It can be seen that results from 
the two approaches are very close to each other. 

3.4.1 Contour Plots of Variables in the X-Y Center Plane 
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Figure 36. — Contours of axial velocity and O 2 mass fraction of the dye in X-Y center plane. 

3.4.2 Radial Profiles of Concentration and Velocity Components at Downstream Locations 


RANS-NL model with & without scalar APDF 
X = 51 mm 




Figure 37. — Radial distributions of concentration at downstream locations. 
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RANS-NL model with & without scalar APDF 
X = 51 mm 



RANS-NL model with & without scalar APDF 



RANS-NL model with & without scalar APDF 




Figure 38. — Radial distributions of velocity components U, Wat downstream locations. 


3.5 Results of URANS Simulation Invoking the Scalar APDF Equation 

In this section, the results of an unsteady RANS (URANS) simulation invoking the scalar APDF 
equation are presented. The scalar APDF equation is solved by an Eulerian method. The main focus is to 
compare the performance of this hybrid approach with the performance of the pure URANS approach. 

Figure 39 shows the contour plots of the mean axial velocity and the mean 0 2 mass fraction of the 
dye in the X-Y center plane, and they are almost identical to their URANS counterparts presented in 
Section 3.2. Figure 40, Figure 41 and Figure 42 present the comparisons of the results from the hybrid 
URANS/APDF simulation and from the pure URANS simulation along the centerline as well as at 
downstream locations x = 25 and 5 1 mm. It can be seen that results from the two approaches are very 
close to each other. 
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3.5.1 Contour Plots of Variables in the X-Y Center Plane 
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Figure 39. — Contours of axial velocity and O 2 mass fraction of the dye in X-Y center plane. 


3.5.2 Centerline and Radial Profiles of Concentration 




Figure 40. — Centerline distributions of concentration and axial velocity. 
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Figure 41. — Radial distributions of concentration at downstream locations. 
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3.5.3 Radial Profiles of Velocity Components at Downstream Locations 


Axial mean velocity profile 





Radial mean velocity profile 



Azimuthal mean velocity 




Figure 42. — Radial distributions of velocity components U, V, W at downstream locations. 
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3.6 Results of TFNS Simulation Invoking the Scalar DWFDF Equation 

In this section, the results of a TFNS simulation invoking the scalar DWFDF equation are presented. 
The scalar DWFDF equation is solved by an Eulerian method. The main focus is to compare the 
performance of this hybrid approach with the performance of the pure TFNS approach. 

Figure 43 shows the contour plots of the instantaneous axial velocity and the instantaneous 0 2 mass 
fraction of the dye in the X-Y center plane. Figure 44, Figure 45 and Figure 46 present the comparisons of 
the results from the hybrid TFNS/DWFDF simulation and from the pure TFNS simulation along the 
centerline as well as at several downstream locations. It can be seen that results from the two approaches 
are very close to each other. 

3.6.1 Contour Plots of Variables in the X-Y Center Plane 
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Figure 43. — Contours of instantaneous axial velocity and instantaneous O 2 mass fraction of the dye in X-Y 
center plane. 


3.6.2 Centerline and Radial Profiles of Mean Concentration 


TFNS, dt=2.e-06, 2nd=0., 4th=0.05, cfl=1.0 



TFNS, dt=2.e-06, 2nd=0., 4th=0.05, cfl=1.0 



Figure 44. — Centerline distributions of concentration and axial velocity. 
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TFNS, dt=2.e-06, 2nd=0., 4th=0.05, cfl=1.0 




Figure 45. — Radial distributions of concentration at downstream locations. 

Radial Profiles of Mean Velocity Components at Downstream Locations 
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Figure 46. — Radial distributions of velocity components U, V, Wat downstream locations. 
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TFNS, Azimuthal mean velocity TFNS, Azimuthal mean velocity 



4.0 Conclusion 

Two groups of validations have been performed using the experiments of a confined swirling coaxial 
water jets flow. The first group focuses on the turbulent scalar flux model to explore the performance of 
the linear formulation versus the nonlinear formulation. Simulations conducted for this group include 
RANS, URANS and TFNS. The second group focuses on the hybrid approach to explore the performance 
of the newly introduced APDF and DWFDF equation and the Eulerian solver for this equation. 
Simulations conducted for this group include RANS/APDF, URANS/APDF and TFNS/DWFDF. 

Regarding the scalar flux model, the linear and nonlinear model have the same behave in RANS and 
URANS simulation. Their respective results shown in various contour plots and profiles at different 
downstream locations are almost identical and in reasonable agreement with the experimental data. In the 
case of TFNS simulations, the differences of results from the linear and nonlinear models are small but 
noticeable. Furthermore, the TFNS results demonstrate significant improvements over their RANS and 
URANS counter parts, when compared with the experimental data. 

Regarding the hybrid approach, the results of RANS/APDF, URANS/APDF and TFNS/DWFDF 
simulations show that they are very close to their respective RANS, URANS and TFNS counterpart. Our 
experiences so far have been that the hybrid approach appears to be more robust in the unsteady 
calculations and to converge faster, hence, using less computing time. 
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Appendix. — Simulations Using Standard k - e Stress 
Model Versus Using Nonlinear Stress Model 

Here we present a set of RANS simulations for the purpose of comparing the performance of the 
standard k - 8 model versus the performance of the nonlinear model. 

1 . RANS simulation with standard k-s model: 

a. Linear algebraic model for Reynolds stresses. 

b. Linear scalar (i.e., heat and species) flux models. 

c. Standard k - 8 transport equations. 

d. Constant = 0.09. 

e. Standard wall function. 

2. RANS simulation with nonlinear stress model: 

a. Nonlinear (up to cubic powers) algebraic model for Reynolds stresses. 

b. Nonlinear scalar (i.e., heat and species) flux models. 

c. Standard k - 8 transport equations. 

d. Variable C ^ based on the realizability constraints. 

e. Generalized wall function including the effects of wall pressure gradients. 

3. Simulations start with the same initial condition, i.e., a RANS solution at the 120,000 th iteration step, 
and use the same numerical parameters: cfl = 1.0, 2nd= -0.01, 4th = 0.05. Both simulations require 
about 40,000 iterations to reach the asymptotic steady state. 

In the following, the presented results include the contour plots of axial velocity U 9 0 2 mass fraction 
of the dye, turbulent kinetic energy k and effective viscosity p + p r , and also the centerline distribution of 
axial velocity U and concentration of the dye, the radial profiles of concentration at the downstream 
locations x = 13, 25, 51 and 102 mm, and the radial profiles of axial velocity U at the downstream 
locations x = 51 and 102 mm. It is quite clear that the standard k- 8 model is not appropriate for 
calculating swirling coaxial jet flow, while the nonlinear stress model performs reasonably well for this 
type of flows. 


RANS standard k-eps model with WF 
cfl=1.0, 2nd=-.01 , 4th=.05, conv=1.e-03 
converged at 155,224 iteration 



u 

■ 3 . 100E+01 
2 . 711E+01 
— I 2 . 322E + 01 
■ 1 . 933E+01 

m 1 . 544E+01 
1.156E+01 
7 . 667E+00 
3 . 778E+00 
I J-l. 111E-01 
™-4. OOOE+OO 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
RANS-NL-wf, cfl=1.0, 2nd=-.01 , 4th=. 05 
At 153,224 iteration, converged Y 



x 


■ 


u 

3 . 100E+01 
2 . 711E+01 
2 . 322E+01 
1 . 933E+01 
1 . 544E + 01 
1 . 156E + 0 1 
7 . 667E+00 
3 . 778E+00 
-1 . 111E-0 1 
-4 . 000E+00 
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RANS standard k-eps model with WF 
cfl=1 .0, 2nd=-.01 , 4th=.05, conv=1.e-03 
converged at 155,224 iteration 


Y 



X 


02 red 

2 . 300E-01 
2 . 044E-01 
1 . 789E-01 
1 . 533E-01 
1 . 278E-01 
1 . 022E-01 
7 . 667E-02 
5 . 111E-02 
2.556E-02 
0 . 000E+00 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 Y 

RANS-NL-wf, cfl=1.0, 2nd=-.01 , 4th=.05 
At 153,224 iteration, converged 


02 red 

h 2 . 30 0 E- 01 
I 2 . 044E-01 
m 1.789E-01 
| 1 . 533E-01 

I 1 . 278E-01 
I 1 . 022E-01 

U 7 . 667E-02 
5 . 11 IE- 02 
2 . 556E-02 
0 . 000E+00 



turb_k 

mm 4 . 053E+01 
I 3 . 603E+01 
■1 3.152E+01 
I 2 . 702E+01 
2 . 252E+01 
I 1.801E+01 

U 1.351E+01 
9 . 007E+00 
4 . 503E+00 
1 . 000E-04 


RANS standard k-eps model with WF 
cfl=1 .0, 2nd=-.01 , 4th=.05, conv=1.e-03 Y 
converged at 155,224 iteration 


turb_k 

4 . 053E+01 
3 . 603E+01 
3.152E+01 
2 . 702E+01 
2 . 252E+01 
1 . 801E+01 
1 . 351E+01 
9 . 007E+00 
4 . 503E+00 
1 . 000E-04 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 Y 

RANS-NL-wf, cfl=1.0, 2nd=-.01, 4th=.05 
At 153,224 iteration, converged 


1 




RANS standard k-eps model with WF 
cfl=1.0, 2nd=-.01 , 4th=.05, conv=1.e-03 y 
converged at 155,224 iteration 


mu 

1.729E-02 
1 . 537E-02 
1 . 345E-02 
1 . 153E-02 
9 . 614E-03 
7 . 6 95E-03 
5 . 776E-03 
3 . 857E-03 
1 . 938E-03 
1.866E-05 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
RANS-NL-wf, cfl=1.0, 2nd=-.01 , 4th=.05 
At 153,224 iteration, converged 



x 


mu 

1 . 729E-02 
1 . 537E- 02 
1. 345E-02 
1 . 153E-02 
9 . 614E-03 
7 . 695E-03 
5 . 776E-03 
3 . 857E-03 
1 . 938E-03 
1 . 866E-05 



RANS, R-J coaxial jet flow 
1 . 5 1- Standard k-eps vs. non-linear model 


O Exp. 

RANS-std-k-eps 
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-100 0 100 200 300 400 500 
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Concentration Concentration 


RANS, R-J coaxial jet flow 




RANS, R-J coaxial jet flow 




RANS, R-J coaxial jet flow 
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